Draft version September 22, 2010 

Preprint typeset using LSTgX style emulateapj v. 1 1/10/09 



A CHANDRA PERSPECTIVE ON GALAXY- WIDE X-RAY BINARY EMISSION AND ITS CORRELATION WITH STAR 
FORMATION RATE AND STELLAR MASS: NEW RESULTS FROM LUMINOUS INFRARED GALAXIES. 

B. D. Lehmer, 1 - 2 D. M. Alexander, 3 F. E. Bauer, 4 - 5 W. N. Brandt, 6 - 7 A. D. Goulding, 3 L. P. Jenkins, 1 - 2 A. Ptak, 2 & 

T. P. Roberts 3 

Draft version September 22, 2010 

ABSTRACT 

We present new Chandra observations that complete a sample of seventeen (17) luminous infrared galaxies 
(LIRGs) with D < 60 Mpc and low Galactic column densities of A/ H ~ 5 x 10 20 cm" 2 . The LIRGs in our 
sample have total infrared (8-1000/im) luminosities in the range of Lir w (1-8) x 10 11 Lq. The high-resolution 
imaging and X-ray spectral information from our Chandra observations allow us to measure separately X-ray 
contributions from active galactic nuclei (AGNs) and normal galaxy processes (e.g., X-ray binaries and hot gas). 
We utilized total infrared plus UV luminosities to estimate star-formation rates (SFRs) and /f-band luminosities 
and optical colors to estimate stellar masses (M t ) for the sample. Under the assumption that the galaxy-wide 
2-10 keV luminosity (Z|£) tr aces the combined emission from high mass X-ray binaries (HMXBs) and low 
mass X-ray binaries (LMXBs), and that the power output from these components are linearly correlated with 
SFR and M*, respectively, we constrain the relation = aM t + /3SFR. To achieve this, we construct a 
Chandra-based data set composed of our new LIRG sample combined with additional samples of less actively 
star-forming normal galaxies and more powerful LIRGs and ultraluminous infrared galaxies (ULIRGs) from 
the literature. Using these data, we measure best-fit values of a = (9.05 ±0.37) x 10 28 ergs s" 1 Mq 1 and 
f3 = (1.62 ±0.22) x 10 39 ergs s" 1 (M Q yr" 1 ) . This scaling provides a more physically meaningful estimate 
of L^, with «0. 1-0.2 dex less scatter, than a direct linear scaling with SFR. Our results suggest that HMXBs 
dominate the galaxy-wide X-ray emission for galaxies with SFR/M* £5.9 x 10~ n yr" 1 , a factor of «2.9 
times lower than previous estimates. We find that several of the most powerful LIRGs and ULIRGs, with 
SFR/M* £ 10~ 9 yr" 1 , appear to be X-ray underluminous with respect to our best-fit relation. We argue that 
these galaxies are likely to contain X-ray binaries residing in compact star-forming regions that are buried 
under thick galactic columns large enough to attenuate emission in the 2-10 keV band (Nh £ 10 23 cm" 2 ). 

Subject headings: galaxies: starburst — infrared: galaxies — X-rays: binaries — X-rays: galaxies — cosmol- 
ogy: observations 



1. INTRODUCTION 

The most energetic nearby star-forming galaxies are lumi- 
nous infrared galaxies (LIRGs), which are classified as having 
total infrared (8-1000/im) luminosities greater than 10 11 Lq 
(where L Q = 3.9 x 10 33 ergs s" 1 is the bolometric luminosity 
of the Sun). The total infrared emission from most LIRGs 
comprises £ 50-90% of their galactic bolometric power out- 
put and is thought to arise primarily from dust reprocessing 
of obscured UV light that originates from underlying star- 
forming regions (see Sanders & Mirabel 1996 for a review). 
Therefore, the total LIRG infrared power provides an effective 
measure of the galactic star-formation rate (SFR; e.g., Kenni- 
cutt 1998) and has been used extensively to select galaxies 
with active star-formation (e.g., Sanders etal. 2003). 
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In normal star-forming galaxies (i.e., those that are not 
dominated by luminous active galactic nuclei [AGNs]), X-ray 
emission originates from X-ray binaries, supernovae, super- 
nova remnants, hot («0.2-l keV) interstellar gas, and O-stars 
(see, e.g., Fabbiano 1989, 2006 for reviews). Several investi- 
gations have now revealed the presence of a strong correlation 
between the galaxy-wide X-ray emission and total SFR (here- 
after, the X-ray/SFR correlation; see, e.g., Bauer etal. 2002; 
Grimm et al. 2003; Ranalli et al. 2003; Gilfanov et al. 2004a,b; 
Persic etal. 2004; Persic & Rephaeli 2007; Lehmer etal. 
2008). The X-ray/SFR correlation is thought to be driven pri- 
marily by X-ray binary emission, and at high X-ray energies 
(2-10 keV), where the emission intensity from hot interstellar 
gas and young stars decreases sharply, X-ray binaries signif- 
icantly dominate the galaxy-wide X-ray power and therefore 
correlate strongly with SFR (e.g., Persic & Rephaeli 2002). 

Detailed Chandra studies of nearby star-forming late- 
type galaxies and passive early-type galaxies have shown 
that the X-ray point source emission from relatively young 
( £ 100 Myr) high-mass X-ray binaries (HMXBs) and older 
( £ 1 Gyr) low-mass X-ray binaries (LMXBs) correlate well 
with galaxy SFR and stellar mass (MJ, respectively (e.g., 
Grimm etal. 2003; Colbert etal. 2004; Gilfanov 2004). Col- 
bert etal. (2004; hereafter C04) utilized C/ianiira-resolved 
point source populations in 32 nearby spiral and elliptical 
galaxies to measure the galaxy-wide X-ray power due to 
point sources Lxp and characterize its correlation with both 
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TABLE 1 

LIRG Sample List and Basic Properties 
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NOTE. — Basic properties of our complete sample of 17 LIRGs with D < 60 Mpc and rV H $ 5 X 10 20 cm -2 , ordered by distance (given in Col. [7]). Col.(l): Common source name. Col.(2) and (3): Right ascension and declination, respectively. Units 
of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminules, and arcseconds. Col. (4): Semi-major axis in arcminutes. Col. (5): Semi-minor axis in areminutes. Col. (6): Position angle in degrees, measured east from 
north. Values given in Col.(2)-(6) were provided by 2MASS with the exception of NGC 3221, which was given by RC3. Col.(7): Luminosity distance in Mpc, as used by Sanders etal. (2003). Col.(8): Galactic column density in units of 10 20 cm -2 as given 
by Dickey & Lockman (1990). Col.(9): Logarithm of the total infrared (8-1000 pirn) luminosity as given by Sanders etal. (2003). Col.(10): Star-formation rate in units of Mq yr -1 . Col.(ll): Logarithm of the stellar mass in units of Mq. Col.(10) and 
(11) were computed following the methods presented in §2.2. Col.(12): Optical morphological classification. When available, these classifications are primarily based onRC3, with the exception of Arp 299, which was based on the Revised Shapley-Ames 
Catalog Of Bright Galaxies (Sandagc & Tammann 1981). Col.(13): Optica! spectroscopic classification of source type. Col.(14): Literature reference for column (13); 1 - Veron-Ccttv & Veron (1986), 2 - Veilleux et al. (1995), 3 - Ho etal. (1997), 4 = 
Kewley etal. (2001), 5 = Yuan el al. (2010), 6 = Armus etal. (1989), 7 =Baan etal. (1998), 8 = Coziol el al. (1998), 9 = Garcia-Marin etal. (2006) and 10 = Keel etal. (1985). 
" Values taken from RC3. 



SFR and M*. The C04 study, which focused on represen- 
tative nearby galaxies with SFR ss 0.01-10 M yr" 1 and 
M* « 10 8 -10 n M Q , showed that L X p is linearly correlated 
with both SFR and M* and the relative contributions from 
HMXBs and LMXBs are directly dependent on specific-SFR 
(i.e., SFR/M*; a measure of the stellar-growth rate). From this 
perspective, the more commonly used X-ray/SFR correlation 
is likely to have significant scatter due to variations in M* for 
a given SFR. 

A recent study by Persic & Rephaeli (2007; hereafter, 
PR07) noted that ultraluminous infrared galaxies (ULIRGs), 
which have SFR > 100 M Q yr _1 , have 2-10 keV luminosity- 
to-SFR ratios that are ~5 times lower than more typical 
nearby galaxies that have SFR <C 100 yr -1 . PR07 suggest 
that the 2-10 keV emission in ULIRGs is likely to be domi- 
nated by HMXBs with negligible contributions from LMXBs, 
which likely play a more significant role in galaxies with 
SFR < 100 Mq yr" 1 . This trend has also been substanti- 
ated by recent Chandra observations of powerful LIRGs and 
ULIRGs from the Great Observatories All-sky LIRG Survey 
(GOALS; Iwasawa etal. 2009; hereafter, Iw09). Despite the 
improved understanding of how the galaxy-wide 2-10 keV 
emission depends on SFR and M* from C04 and sugges- 
tions that extremely powerful LIRGs and ULIRGs have X-ray 
power dominated by young HMXB populations (from e.g., 
PR07 and Iw09), researchers have yet to test directly and con- 
strain the C04 relation using observations of galaxies with 
SFR > 10 Mq yr" 1 . 

In this paper, we present new Chandra results for a volume- 
limited sample (D < 60 Mpc) of 17 LIRGs (including 10 new 
LIRGs from this study; see Tables 1 and 2) with Galactic col- 
umn densities tVh ~5 x 10 20 cm" 2 and L IR ssl-8 xl0 n L e 
(SFR > 10M Q yr" 1 ; see Figs. 1 and 2). We first utilize the 
high-resolution Chandra data to separate spatially and mea- 
sure contributions from normal galaxy-wide X-ray emission 
(i.e., not due to an AGN) and AGN when present; this enables 



us both to investigate the X-ray/SFR correlation in detail for 
LIRGs and to characterize the AGN activity in the population. 
Using the new LIRG constraints in combination with both the 
Chandra results from C04 for less actively star-forming nor- 
mal galaxies and Iw09 for more actively star-forming LIRGs 
and ULIRGs, we provide improved constraints upon how the 
normal-galaxy 2-10 keV luminosity depends on SFR and M*, 
and therefore the relative contributions from HMXBs and 
LMXBs. This analysis has the advantage over previous stud- 
ies in that it uniquely covers «4 orders of magnitude in SFR 
and ~3.5 orders of magnitude in M* and is uniformly based 
on Chandra data. 

Throughout this paper, we make estimates of SFR and 
M* using a Kroupa (2001) initial mass function (IMF); 
when making comparisons between these estimates and those 
quoted in other studies, we have adjusted all values to cor- 
respond to our adopted IMF. We sometimes make com- 
parisons with the Milky Way (MW). We adopt values of 
M* w 5 x 10 10 M Q (Hammer etal. 2007), SFR « 2.0 M 
yr" 1 (McKee & Williams 1997), and 2-10 keV luminosity 
Lg£ = 3.3 x 10 39 ergs s" 1 (Grimm etal. 2002) for the MW 
H = 70 km s" 1 Mpc" 1 , tt M = 0.3, and Cl\ = 0.7 are adopted 
throughout this paper (e.g., Spergel et al. 2003). 

2. SAMPLE 

2.1. LIRG Sample Selection 

Using the IRAS Revised Bright Galaxies Sample (RBGS; 
Sanders etal. 2003), we constructed a sample of nearby 
LIRGs that were at distances less than 60 Mpc and had Galac- 
tic column densities Nh < 5 x 10 20 cm" 2 (Dickey & Lock- 
man 1990). The former requirement results in Chandra imag- 
ing resolution of ^150 pc in the region of the nucleus, which 
allows for reasonable distinction between nuclear AGNs and 
starburst regions (see, e.g., Levenson et al. 2004 for an anal- 
ysis of the AGN plus starburst nucleus of NGC 5135 in 
our sample). The latter requirement results in moderate-to- 
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FIG. 1 . — Digitized sky survey (DSS) red images of each of the 17 LIRGs that make up our sample. The panel sizes presented here vary to illustrate galaxy 
morphology most clearly, and in the lower left and right corners of each panel, we have provided vertical reference bars with physical sizes of 10 kpc and angular 
sizes of 1 arcmin, respectively. Dashed gray ellipses represent the total A^-band galaxy size and orientation as described in columns (2)-(6) of Table 1 and gray 
crosses indicate the adopted locations for the galactic nuclei (see §4.2 for details). These images show that the galaxies making up our LIRG sample have a 
diversity of optical morphological types, physical sizes, and inclinations. 



good Galactic transparency to soft X-ray emission, which is 
valuable for detecting and modeling the emission from hot 
~0.5-l keV gas that is often found in star-forming galaxies. 
In total, 36 RBGS LIRGs reside at D < 60 Mpc and 17 satisfy 
our Galactic column density criterion. In Table 1, we provide 
the basic properties of our LIRG sample, and in Figure 1, we 
show Digitized Sky Survey (DSS) red images of the galax- 
ies in our sample. In contrast to ULIRGs, which are nearly 
all late-type galaxies with evidence for mergers/interactions 
(e.g., Melnick & Mirabel 1990; Sanders etal. 1988; Chen 
etal. 2010), the LIRGs that make up our sample consist of 
a variety of optical morphological types (e.g., spirals, barred 
spirals, SOs) with only a small fraction (~10%) showing evi- 
dence for mergers/interactions. 

In Figure 2, we show the infrared luminosity Lir versus dis- 
tance for the RBGS and highlight our sample (filled black cir- 
cles). Of the 17 LIRGs in our sample, seven (NGC 1068, 
NGC 1365, NGC 7552, NGC 4418, NGC 4194, Arp 299, 
and NGC 5135) had Chandra archival data. Detailed Chan- 
dra analyses for these sources can be found in the literature 
for NGC 1068 (Young etal. 2001; Smith & Wilson 2001), 
NGC 1365 (Risaliti etal. 2009a,b; Soria etal. 2009; Strat- 
eva & Komossa 2009; Wang et al. 2009), NGC 4194 (Kaaret 
& Alonso-Herrero 2008), Arp 299 (Zezas etal. 2003), and 
NGC 5135 (Levenson etal. 2004). To enable new Chan- 
dra constraints on the LIRG population in general, we pro- 
posed successfully in Chandra Cycle 10 (PI: D. M. Alexan- 
der) to observe the remaining 10 LIRGs. Our Chandra pro- 
gram (and consequently adopted exposure times; see Ta- 
ble 2) was designed to allow for the detection of ultralumi- 
nous X-ray sources (ULXs) down to a 0.5-8 keV limit of 



«2 x 10 39 ergs s" 1 for the 17 LIRGs. In a forthcoming paper 
(Jackson et al. in-prep), we will present the X-ray and statis- 
tical properties of the ULX populations detected in the LIRG 
sample. 

As is evident from Figure 2, our sample provides a use- 
ful complement to the Lir distribution of Chandra-observed 
galaxies between less active nearby galaxies (e.g., from C04) 
and more powerful LIRGs and ULIRGs (e.g., from the 
GOALS; Iw09). We note that the RBGS sources with Chan- 
dra observations highlighted in Figure 2 do not constitute a 
complete census of the Chandra archive. However, incom- 
pleteness is only most severe for galaxies with Lir £ 10 n L Q , 
which are not the primary focus of this analysis; detailed 
Chandra and multiwavelength studies of such galaxies will 
be the subject of forthcoming work by the Spitzer Infrared 
Nearby Galaxies Survey (SINGS) collaboration (Jenkins et al. 
in-prep; Ptak etal. in-prep; see Kennicutt etal. 2003 for sur- 
vey details). For LIRGs and ULIRGs, the Chandra observa- 
tions highlighted in Figure 2 are nearly complete. 



2.2. Stellar Mass and Star Formation Rate Estimates 

We estimated stellar masses for the LIRGs in our sample 
following the prescription outlined in Appendix 2 of Bell et al. 
(2003), which provides galaxy mass-to-light ratios in several 
optical/near-IR bands as a function of optical/near-IR color. 
We utilized B - V colors from the NAS A/IPAC Extragalactic 
Database (NED), which were primarily from the Third Refer- 
ence Catalog of Bright Galaxies (RC3; de Vaucouleurs etal. 
1991). TT-band luminosities were computed using data from 
the Two Micron All Sky Survey (2MASS) large galaxy atlas 
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FIG. 2. — Total infrared (8-1000/^m) luminosity Lir versus distance for 
IRAS sources (light-gray filled circles), as provided by Sanders etal. (2003). 
Our complete sample of CAanrfra-observed LIRGs that have D < 60 Mpc 
and Nu £ 5 X 10 10 cm -2 (dashed boundary lines) are indicated with black 
filled circles. Nearby normal late-type galaxies that have been studied ex- 
tensively with Chandra have been highlighted with open triangles (C04) 
and LIRGs/ULIRGs from the Great Observatories All-sky LIRG Survey 
(GOALS) sample (Iw09) are indicated with open squares. Our sample pro- 
vides new constraints in the range between nearby normal galaxies and pow- 
erful LIRGs/ULIRGs. 

(Jarrett etal. 2003) and ZQ-band extended source catalogs^ 
We note that for NGC 3221, no 2MASS counterpart was avail- 
able, and we therefore adopted the estimate provided by 
the GOALS team (Howell etal. 2010; hereafter, H10), ad- 
justed to match our adopted distance and IMF (see below). 
For the 16 LIRGs having 2MASS counterparts, we computed 
stellar masses using the following equation: 

logM*/M Q = logL jC /L^. o + 0.135(B-y)-0.356. (1) 

The numerical constants in equation 1 were supplied by Ta- 
ble 7 of Bell et al. (2003) and are appropriate for our choice 
of B— V color and Lr-; the normalization has been adjusted by 
0.15 dex to account for our adopted Kroupa (2001) IMF. Un- 
certainties at the «0. 1 dex level are expected due to variations 
in stellar ages, dust attenuation, and bursts of star formation 
(see Bell et al. 2003 for further details). 

We note that 12 of the LIRGs in our sample (including 11 
LIRGs with stellar mass estimates from equation 1) have been 
studied extensively from the far-UV to IR (including observa- 
tions with both GALEX and Spitzer) by the GOALS collabo- 
ration and have stellar mass estimates available (e.g., Armus 
etal. 2009; H10). After accounting for known differences in 
adopted distances and IMF, we find that our stellar mass es- 
timates are in good agreement with those of H10, with our 
values being larger by w4% (median value) with a la scatter 
of«0.09dex. 

By selection, our LIRG sample consists of galaxies radi- 
ating powerfully over the far-IR (8-1000 /im) wavelength 
range. The far-IR light is expected to be produced by dust- 
reprocessed emission from underlying UV-obscured star- 
formation activity and can therefore be used as a direct mea- 
surement of the SFR. Using the 8-1000 /im total infrared lu- 
minosities Lir provided by Sanders et al. (2003) and the meth- 
ods described in § 3.2 of Bell etal. (2005), we estimated the 

8 The 2MASS extended source catalog is available online via 
http://tdc-www.harvard.edu/catalogs/tmxsc.html 



total galaxy SFR using the following equation: 

SFR(M vr" 1 ) = 7 9.8 x lO" 11 ^, (2) 

where Lir is expressed in units of solar bolometric luminosity 
(Lq = 3.9 x 10 33 ergs s" 1 ) and the factor 7 = 1.000-1. 024 pro- 
vides small corrections that account for UV emission emerg- 
ing from unobscured star-forming regions. Equation 2 is a 
direct variant of equation 1 presented by Bell et al. (2005) and 
was derived from PEGASE stellar-population models, which 
assumed a 100 Myr old population with constant SFR and a 
Kroupa (2001) IMF. This calibration uses identical assump- 
tions to the SFR calibrations provided by Kennicutt (1998). 
For the 12 LIRGs that overlap with the H10 sample, 7 was 
provided by H10 using values of the total infrared flux from 
IRAS (/ir) and the far-UV flux (/prjv) from GALEX: 7 = 
1+/fuv//ir- For these 12LIRGs, we find good agreement be- 
tween SFR measurements, with our SFR values being larger 
by «4% (median offset) with small scatter («0.003 dex; la). 
For the five LIRGs in our sample without H10 counterparts, 
we adopted ^ medmn « 1.017, the median value of 7 for the 12 
LIRGs with H10 counterparts. 

In Table 1, we provide the derived stellar masses and SFRs 
for our LIRGs. The sample spans w (1.3-16) xl0 1() M Q 
(median value of M* « 7.2 x 1O 1O M ) and SFR = 9.8- 
75.6 M Q yr" 1 (median value of SFR « 11 M Q yr" 1 ). For 
comparison, the LIRGs/ULIRGs from the Chandra GOALS 
project (Iw09) span a similar range of stellar masses (M* ~ 
(1.6-27) xl0 10 M Q ; median value of M* sa 6.5 x 1O 1O M ), 
but contain more actively star-forming galaxies (SFR w 
14-215 M yr" 1 ; median value of SFR w 68 M yr" 1 ). The 
nearby late-type galaxies in the C04 study span the stellar 
mass range of «6 x 1O 7 -1O H M (median value of M+ w 
1.8 x 1O 1O M ) and less powerful star-forming galaxies with 
SFR = 0.03-5.6 M Q yr" 1 (median value of SFR = 0.8 M 

yr" 1 ). 

3. CHANDRA OBSERVATIONS AND DATA REDUCTION 

In Table 2, we present the Chandra ACIS observation log 
for our LIRG sample. Our analyses began with the Chandra 
X-ray Center (CXC) processed Level 2 events files, which 
were all processed using pipeline version 7.6.8 or greater. 
Additional reductions and analysis of the data were per- 
formed using the Chandra Interactive Analysis of Observa- 
tions (CI AO) Version 4.1 tools and custom software. We 
screened our events files for undesirable grades using the stan- 
dard ASCA grade set (ASCA grades 0, 2, 3, 4, 6) and ex- 
cluded from the events lists background flaring events that 
were ^ 3a times higher than normal. For the observations 
344 (NGC 1068), 3554 (NGC 1365), 7071 (NGC 4194) and 
2187 (NGC 5135) we removed short flaring events lasting 
«0.5, 1.0, 0.5, and 3.0 ks, respectively. 

Out of the 17 LIRGs in our sample, three galaxies 
(NGC 1068, 1365, and 4418) have more than one Chandra 
observation available. However, in the case of NGC 1365, 
observations 6868, 6869, 6870, 6871, 6872, and 6873 cover 
only a small fraction of the entire galactic extent as defined 
in Table 1 and constitute ss86.5% of the total ACIS exposure 
available. In order to avoid galaxy-wide spectral fits (see § 4. 1 
below) that are statistically biased towards these smaller non- 
representative regions, we therefore chose for this study to 
make use of only observation 3554, which covers the entire 
galactic extent and is sufficiently sensitive for our study. 
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TABLE 2 
Chandra OBSERVATION LOG 



Source Name 


OBS. ID 


Obs. Start 
(UT) 


Exposure Time" 
(ks) 


Obs. Mode* 


Camera 


SUBARRAY? C 
(Y/N) 


Pipeline Version^ 


NGC 1068 


344 


2000 Feb 21, 15:48 


46.9 


F 


ACIS-S 


N 


7.6.11.2 




370 


2000 Feb 22,05:51 


3.7 


VF 


ACIS-S 


Y 


7.6.11.1 


NGC 1365 


3554 


2002 Dec 24, 14:57 


13.6 


F 


ACIS-S 


N 


7.6.9 




6868 c 


2006 Apr 17, 19:05 


14.6 


F 


ACIS-S 


Y 


7.6.7.1 




6869 e 


2006 Apr 20, 10:00 


15.5 


F 


ACIS-S 


Y 


7.6.7.1 




6870 c 


2006 Apr 23, 10:05 


14.4 


F 


ACIS-S 


Y 


7.6.7.2 




6871 c 


2006 Apr 10, 07:02 


13.2 


F 


ACIS-S 


Y 


7.6.7.1 




6872 e 


2006 Apr 12, 12:31 


14.6 


F 


ACIS-S 


Y 


7.6.7.1 




6873 e 


2006 Apr 14, 23:26 


14.6 


F 


ACIS-S 


Y 


7.6.7.1 


NGC 7552 


7848 


2007 Mar 31, 05:48 


5.1 


F 


ACIS-S 


N 


7.6.11 


NGC 4418 


4060 


2003 Mar 10, 00:58 


19.8 


F 


ACIS-S 


Y 


7.6.8 




10391 


2009 Feb 20, 20:34 


5.7 


VF 


ACIS-S 


Y 


7.6.11.10 


NGC 4194 


7071 


2006 Sep 09, 23:30 


35.0 


F 


ACIS-S 


N 


7.6.8.1 


IC 5179 


10392 


2009 Jun21, 13:51 


12.0 


VF 


ACIS-S 


Y 


7.6.11.10 


ESO420-G013 


10393 


2009 May 13, 05:46 


12.4 


VF 


ACIS-S 


Y 


7.6.11.10 


Arp 299 


1641 


2001 Jul 13, 11:07 


24.3 


F 


ACIS-I 


Y 


7.6.9 


NGC 838 


10394 


2008 Nov 23, 09:03 


13.8 


VF 


ACIS-S 


Y 


7.6.11.9 


NGC 5135 


2187 


2001 Sep 04, 15:23 


26.4 


F 


ACIS-S 


Y 


7.6.9 


NGC 5394/5 


10395 


2009 Mar 18, 23:28 


15.7 


VF 


ACIS-S 


Y 


7.6.11.10 


NGC 5653 


10396 


2009 Apr 11, 05:46 


16.5 


VF 


ACIS-S 


Y 


7.6.11.10 


NGC 7771 


10397 


2009 May 22, 05:12 


16.7 


VF 


ACIS-S 


Y 


7.6.11.10 


NGC 3221 


10398 


2009 Mar 19,04:19 


19.0 


VF 


ACIS-S 


Y 


7.6.11.10 


CGCG 049-057 


10399 


2009 Apr 17, 05:43 


19.1 


VF 


ACIS-S 


Y 


7.6.11.10 


IC 860 


10400 


2009 Mar 24,02:38 


19.1 


VF 


ACIS-S 


Y 


7.6.11.10 


NGC 23 


10401 


2008 Oct 27, 06:22 


19.4 


VF 


ACIS-S 


Y 


7.6.11.9 



NOTE. — Links to the data sets in this table have been provided in the electronic edition. 

" All observations were continuous. These times have been corrected for removed data that was affected by high background; see § 3. 
* The observing mode (F=Faint mode; VF=Very Faint mode). 
c States whether a subarray was used in the observation. 

d The version of the CXC pipeline software used for basic processing of the data. 

e Observation covers only a small fraction of the galactic extent and was therefore not used in the analyses presented in this paper. 



For the remaining two LIRGs with more than one Chan- 
dra observation (NGC 1068 and 4418), we corrected the rel- 
ative astrometry to match the frames with the longest expo- 
sures. We first ran wavelet ect at a false-positive proba- 
bility threshold of 10" 6 to create point-source catalogs for 
each frame. For a given galaxy, we filtered each of our 
wavdetect catalogs to include only sources that were com- 
mon to all frames and within 6' of the Chandra aim-point, 
where the point-spread function size and resulting positional 
errors are small. Using these source lists, we registered each 
aspect solution and events list to the frame with the longest 
exposure time using the CIAO tools repro ject_aspect 
and repro ject_events, respectively. The resulting as- 
trometric reprojections gave nearly negligible linear trans- 
lations (<0.38 pixels), rotations (<0.1 deg), and stretches 
(<0.06% of the pixel size). Using the astrometrically- 
reprojected events lists, we combined the observations using 
dmmerge to create merged events lists. 

Using the events lists, we constructed images of each 
LIRG in three bands: 0.5-8 keV (full band; FB), 0.5-2 keV 
(soft band; SB), and 2-8 keV (hard band; HB). For each 
of the three bands, we constructed corresponding exposure 
maps following the basic procedure outlined in §3.2 of Horn- 
schemeier etal. (2001); these maps were normalized to the 
effective exposures of sources located at the aim points. 
This procedure takes into account the effects of vignetting, 
gaps between the CCDs, bad column and pixel filtering, the 
spatially-dependent degredation of the ACIS optical blocking 



filter, and the reduced area of frames that use subarrays (see 
Table 2). A photon index of Y = 1.4, the slope of the ex- 
tragalactic cosmic X-ray background in the full band (e.g., 
Hickox & Markevitch 2006), was assumed in creating the ex- 
posure maps. 

4. ANALYSES 

4.1. Extraction of Total LIRG X-ray Properties 

To measure the total galaxy-wide X-ray emission originat- 
ing from each of the LIRGs, we began by extracting total 

galactic counts Sf^ and effective exposure times T^f from 
each of the three images and exposure maps using elliptical 
regions that approximate the galactic stellar extent. When- 
ever possible, we adopted for these elliptical regions the "to- 
tal" semi-major/minor axes and position angles provided by 
the Two Micron All Sky Survey (2MASS) large galaxy at- 
las (Jarrett et al. 2003) and the 7sT s -band extended source cat- 
alogs; these regions, which were determined using a curve- 
of-growth analysis (see §4.2 of Jarrett etal. 2003), are shown 
in Figure 1 (dashed ellipses) and their values are tabulated 
in Table 1. For NGC 3221, 2MASS values were not avail- 
able and we therefore adopted the RC3 values for the semi- 
major/minor axes and position angle for this galaxy. We note 
that for 14 of the 17 LIRGs in our sample, both 2MASS and 
RC3 galaxy extent parameters were available. We find that 
the 2MASS values of major and minor axes are on average 
larger than those of RC3 by only small factors of «1.13 and 
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TABLE 3 

Basic X-ray Properties and Spectral Fitting Results 





Net Counts 














Other 


AGN 










kT 






log< 




Evidence 


Dominant? 


Source Name 


(0.5-8 keV) 


(2-8 keV) 


r 


(keV) 


X 2 /" 


V 


(ergss ') 


^2-8 keV/ v 2-8 keV 


for AGN 


(Y/N) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 1068 


150875 ± 404 


11215 ± 131 










41.49" 


0.74" 


Opt., IR, X-ray 


Y 


NGC 1365 


4719 ± 127 


841 ± 90 


1.88 ±0.12 


0.73 ± 0.05 


1.26 


84 


40.65 


0.88 


Opt., IR, X-ray 


Y 


NGC 7552 


939 ± 40 


113 ± 23 


2.06 ±0.13 


0.64 ± 0.04 


0.91 


46 


40.40 


<0.08 




N 


NGC 441 8 


49 ±23 


<53 


3.28 ± 0.46 








<39.58 






N 


NGC 4194 


2413 ± 53 


263 ± 24 


2.11 ±0.08 


0.83 ± 0.02 


1.23 


79 


40.43 


0.50 


IR 


Y 


IC5179 


520 ± 33 


62 ± 22 


2.10 ±0.19 


0.84 ± 0.06 


0.62 


32 


40.49 


<0.09 




N 


ESO420-G013... 


762 ±31 


53 ± 15 


3.00 ±0.19 


0.71 ± 1.61 


1.09 


30 


40.08 


0.43 


Opt., IR 


N 


Arp 299 


4456 ± 72 


836 ± 36 


1.85 ±0.06 


0.80 ± 0.02 


1.29 


130 


41.32 


0.11 


Opt., IR, X-ray 


N 


NGC 838 


632 ± 29 


70 ± 14 


2.07 ±0.13 


0.73 ± 2.84 


0.83 


28 


40.61 


0.41 




N 


NGC 5135 


3756 ± 74 


346 ± 37 


2.53 ± 0.08 


0.78 ± 0.55 


1.35 


117 


40.70 


0.44 


Opt., IR, X-ray 


N 


NGC 5394/5 


289 ± 38 


92 ± 30 


1.69 ±0.38 


0.09 ± 0.20 


0.65 


37 


40.49 


0.27 


Opt. 


N 


NGC 5653 


389 ± 30 


50 ± 19 


2.21 ± 0.26 


0.68 ± 0.06 


1.20 


26 


40.18 


<0.11 




N 


NGC 7771 


864 ± 42 


233 ± 30 


1.24 ±0.21 


0.63 ± 2.39 


1.03 


55 


41.18 


0.53 


IR, X-ray 


Y 


NGC 3221 


310 ±31 


69 ±22 


1.58 ±0.16 


fs0.0 r 


0.95 


26 


40.66 


0.40 




N 


CGCG 049-057. .. 


33 ± 13 


<29 


1.05 ±0.57 








<40.29 






N 


IC 860 


<44 


<35 


R*2.0* 








<40.19 






N 


NGC 23 


717 ±38 


<71 


2.03 ± 0.21 


0.73 ± 1.12 


1.18 


41 


<40.54 


>0.24 




N 



NOTE. — Basic X-ray properties and spectral fitting results for the LIRGs in our sample. Col.(l): Source name. Col. (2) and (3): Net source counts in the 
0.5-8 keV and 2-8 keV bandpasses, respectively (see §4.1 for details). For sources with ^,300 FB counts listed in column 2, X-ray spectral fitting was 
performed over the 0.5-8 keV band. Col. (4) and (5): Best-fit values of T and kT, respectively, for an assumed power-law plus Raymond-Smith thermal plasma 
model (see §4.1 for details). Col.(6) and (7): Reduced \ 2 value (x 2 M and the number of degrees of freedom in the fits (u), respectively. Col. (8): Integrated 
2-10 keV total luminosity, computed using our best-fit models and distances provided in Table 1, column (7). Col.(9): Ratio of 2-8 keV count-rates from the 
nuclear region of each galaxy and the total galaxy- wide extent. Count-rates were computed following the methods described in §§4.1 and 4.2. A source with 
large nuclear contributions to the total 2-8 keV count-rates harbors either an AGN or a strong nuclear starburst. Col. (10): Notes whether there is evidence for 
AGN activity via optical spectroscopy ("Opt"; see columns 13 and 14 in Table 1), the mid-infrared spectrum ("IR"), or from X-ray studies in the literature 
("X-ray"); see § 4.2. Col.(ll): Indicates whether the integrated 2-10 keV luminosity (column 8) is dominated by AGN activity. We consider sources that have 

both <J>™8 kl; v/*2-8 keV — °- 5 (column 9) and additional evidence for AGN activity from column 10 to have the 2-10 keV luminosity dominated by an AGN 
("Y"); otherwise, we consider the galaxy to be dominated by X-ray binary emission ("N"). 
"Values are inferred using detailed analysis by Young etal. (2001). 
b A simple power-law model with T 2.0 was assumed for IC 860. 
r Best-fit spectrum did not require the use of a hot gas component. 



wl.08, respectively. 

To estimate the background contributions to the total counts 
for each galaxy, we extracted background counts Sf^ g and 

exposure times 7^ for each of the three images using 3-5 
moderate size («20-90 arcsec radius) circular regions that 
were outside the galactic 7T v -band ellipses and did not over- 
lap within 10" of the galaxy or any sources detected by 
wavdetect at a false-positive probability threshold of 10~ 6 . 
Using the galaxy and background counts, we estimated the 

net galaxy-wide LIRG counts JV gal = Sf^-S^ g T^ /T^ for 

each of the three bands. Since T^f and 7j| *J are computed us- 
ing the sum of exposure map values over the source and back- 
ground regions, respectively, the estimated background con- 
tribution to the net counts (i.e., S^ T^f/T^) are corrected 
for both extraction area and vignetting. We note that while 
these background estimates will account for emission due 
to unrelated background X-ray sources (e.g., distant AGNs; 
see, e.g., Morretti etal. 2003; Bauer etal. 2004; Kim etal. 
2007) that fall below our source detection threshold, these es- 
timates do not include contributions from unrelated sources 
that are detected in the X-ray band that may lie within the 
extent of the galaxy. Using the best-fit logAMogS from Mor- 
retti et al. (2003) and the elliptical galaxy areas described in 
columns (2)-(6) of Table 1, we estimate that unrelated X-ray 



detected sources are expected to contribute ^,8% (median 
«2%) to the total 2-8 keV band count-rate of each LIRG in 
our sample. Since these contributions are small and subject to 
large uncertainties, we therefore do not make any additional 
corrections to the net galaxy count-rates. 

We considered a source to be detected in a given band if 
the signal-to-noise ratio (S/N = N s ' di /[S^ g T t tf /T^] 05 ) was 
greater than or equal to three. When a source was not de- 
tected in a given band, we computed 3a upper limits on its net 
counts. In columns (2) and (3) of Table 3, we provide the net 
0.5-8 keV and 2-8 keV source counts for the 17 LIRGs in our 
sample, and in Figure 3, we show their 0.5-8 keV adaptively- 
smoothed images. We note that all galaxies, with the excep- 
tion of IC 860, are detected in the full band. In the hard band, 
NGC 4418, CGCG 049-057, IC 860, and NGC 23 are un- 
detected. The LIRGs have net full-band source counts in the 
range of «30-10 5 counts; for the 16 sources that were de- 
tected in the full band, X-ray spectral fitting was performed. 

To model the basic spectral properties of the full-band de- 
tected LIRGs in our sample, we utilized the X-ray spectral 
modelling program xspec Version 12.5.1 (Arnaud 1996). 
We note that the galaxy-wide and nuclear spectral properties 
of the LIRG NGC 1068 have been studied extensively in the 
literature by, e.g., Young et al. (2001); hereafter, we adopt the 
complex best-fit spectral models and resulting properties pro- 



X-ray Emission from LIRGs 7 



NGCJQ66 

/ s 
\ % I 


NQPJ365 

* . \ 




NGC 7553 

\ I 

."" 1 


NGC 4416 

\ + s 




NGC 4194 

i r * 

A 


IC 5179 

1 ^ " " 1 


ESO 430-G013 


Arp 399 

^ - s 
/ \ 

1 #>' 1 

v - f 

i T 




NGC 838 

• 

• -- — -~ 
I 1 ▼ ' 

N __ _ ' 




NGC 6135 

■ / ' 

s 

l * 




NGC 5394/5 

t \ 

\*\ 

V, * < 1 
1 1 


NGC 5653 ■ 

1 -•' i 


NGC 7771 

*- — ^ 

1. *** 

E 1 < 1 — 


NCC 3331 

I ' ■ 

J*' 

(1 

1 ' 




CGCC 49-05* 

r * 

\ ♦ * 

v » 




IC 860 

f \ 
S + I 

* 




NGC 33 

/ « 

' / 
1 ^ 

* ♦ 





FIG. 3. — Full band (0.5-8 keV) adaptively smoothed images of the 17 LIRGs in our sample. The images were created using the CIAO tool csmooth, applied 
to the raw images at the 2.5(7 level. Symbols and scales have the same meaning as they did in Figure 1; however, the image sizes have been adjusted to show 
most clearly the X-ray morphologies of the galaxies. Galaxies with 2-8 keV emission dominated by an AGN component (see § 4.2) have been denoted with an 
"A" in the lower right. We note that each image has a different exposure time (see Table 2) and the smoofhed-image sharpness generally increases with increasing 
exposure time. 



vided by Young et al. (2001). 

For each observation of the remaining 15 galaxies with full- 
band detections, we employed the CIAO tool specextract 
to construct source and background spectra, response ma- 
trix files (RMFs), and properly-weighted ancillary response 
files (ARFs) over the source and background regions de- 
scribed above. For the 13 full-band detected LIRGs that 
had 300 net counts (excluding NGC 1068), we utilized 
X 2 fitting of the data grouped in bins with a minimum of 
20 counts. For these galaxies, we assumed a simple spec- 
tral model consisting of Galactic absorption (as given in col- 
umn [8] of Table 1), a Raymond-Smith thermal plasma (from 
hot gas associated with star-formation activity) with solar 
abundances, plus a power-law component (originating from 
X-ray binaries and AGN when present). Similar composite 
fits have been performed successfully for ULIRGs (see, e.g., 
Franceschini etal. 2003; Ptak etal. 2003). We find that our 
simple model provides a reasonable fit to these 13 LIRGs 
(X 2 /V= 0.62-1.35; median[x 2 /V] = 1.09). For NGC 3221, 
we find that the best-fit plasma component provides only 
a negligible contribution to the total spectrum; we there- 
fore fit the X-ray spectrum using only a power-law model. 
For the remaining two full-band detected LIRGs that had 
«30-50 net counts (NGC 4418 and CGCG 049-057) we per- 
formed simple power-law fits to the unbinned spectral data by 
minimizing the Cash statistic (Cash 1979). For NGC 4418, 
which was observed in two exposures (see Table 2), we per- 
formed joint-spectral fitting of the two observations. 

The galaxy-wide best-fit model parameters are provided in 
Table 3. Using these models and the distances provided in 
column (7) of Table 1, we computed the total galaxy-wide 



2-10 keV power output Lf^. For sources that were not de- 
tected in the hard band, we computed 3<r upper limits to 
using our best-fit spectral model renormalized to the 3a upper 
limit on the 2-8 keV count-rate. In the case of IC 860, which 
was not detected in the full or hard bands, we assumed a sim- 
ple power-law model with T = 2.0. Our best estimate values 

of L| dl are provided in column (8) of Table 3. Using our spec- 
tral fits, we estimate that hot gas emission will provide only 
a negligible contribution (<6.5% for all LIRGs) to the total 
2-10 keV emission for our LIRGs. We therefore expect that 
LIRGs that do not harbor a luminous AGN (see § 4.2 below) 
will have dominated by X-ray binary emission. 

4.2. AGN Contribution to Total Galaxy X-ray Emission 

As noted in columns (13) and (14) of Table 1, «5-7 of the 
LIRGs in our sample have been identified to harbor an AGN 
via optical spectroscopy. Additional evidence for AGN activ- 
ity in these LIRGs has been observed in Spitzer infrared spec- 
trograph (IRS) data available from the archive for the entire 
sample. These data were analyzed following the methods out- 
lined in Goulding & Alexander (2009). The Spitzer data re- 
veal emission-line features (e.g., [Ne V] A14.3^mor24.3/im) 
indicative of AGN activity for NGC 1068, NGC 1365, 
NGC 4194, ESO 420-G013, NGC 5135, and NGC 7771, and 
in the case of Arp 299 a hot-dust component from an AGN has 
been identified by Alonso-Herrero et al. (2009). Furthermore, 
X-ray studies of NGC 1068 (e.g., Ogle etal. 2003), NGC 1365 
(e.g., Risaliti etal. 2009b), Arp 299 (e.g., Delia Ceca etal. 
2002), NGC 5135 (Levenson etal. 2004), and NGC 7771 
(Jenkins et al. 2005) have revealed X-ray spectral signatures 
of AGN activity (primarily through the presence of a heavily- 
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absorbed X-ray emission component and/or Fe K emission 
lines). In column (10) of Table 3, we highlight whether a 
source has evidence for AGN activity from the optical, in- 
frared, and/or X-ray bands. 

In cases where the AGN is luminous in the X-ray band, the 
AGN component will dominate the integrated (galaxy-wide 
plus nuclear) 2-10 keV power Lf^ (column 8 of Table 3). In 
a forthcoming paper (Stratford et al. in-prep), we will present 
multiwavelength analyses that characterize the detailed ac- 
cretion properties of the AGNs in our LIRG sample; how- 
ever, for our purposes, we are primarily interested in how the 
galaxy-wide 2-10 keV emission from X-ray binaries corre- 
lates with SFR and and therefore remove the sources that 
have 2-10 keV emission dominated by AGN activity. 

To estimate the AGN contribution to for each LIRG, 
we began by estimating the position of the nucleus of each 
galaxy. For the majority of the LIRGs, the nucleus was iso- 
lated first by eye using the DSS red images and then by 
choosing manually the wavdetect position of the bright- 
est 0.5-8 keV source within the optical nuclear region. The 
adopted LIRG nuclei are generally consistent with the appar- 
ent optical nuclei; their locations are shown in Figures 1 and 
3 (gray crosses). In most cases, the X-ray nucleus can be 
obviously characterized as a bright point-source at the center 
of diffuse X-ray emission or a single point-source at the cen- 
ter of the optical nuclear region. However, for Arp 299 and 
NGC 5135, there was more than one possible choice for the 
central nuclear X-ray source; for these galaxies, we adopted 
the positions of the AGN components noted by Zezas etal. 
(2003) and Levenson et al. (2004), respectively. Furthermore, 
for NGC 5653 and IC 860 there were no obvious central X-ray 
point sources to identify as the nucleus; we therefore esti- 
mated by eye the nuclear positions of these galaxies using the 
DSS red images. 

At the nuclear position of each source, we extracted nuclear 
source counts 5™ t c and mean exposure times r t " t uc from the 
2-8 keV images and exposure maps, respectively, using a cir- 
cular aperture with radius equal to the «90% encircled-energy 
fraction radius (Feigelson et al. 2000)0 We estimated the lo- 
cal nuclear background count-rate for each source by extract- 
ing nuclear background counts ,Sj|]j5 and exposure times T b ™ c 
from an annulus centered on the nuclear position with inner 
and outer radii equal to 1.5 and 2.5 times the «90% encircled- 
energy fraction radius, respectively. The net nuclear 2-8 keV 
count-rates $™ c 8 keV = 5™ c /T™ c -S^/T$£ were then com- 
puted and compared with the total galaxy-wide count-rates 
3>2-r keV computed in § 4. 1 . 

In column (9) of Table 3, we provide $™ c 8 keV /$^ 8 keV for 
the LIRGs in our sample. We note that for NGC 1068, we uti- 
lized the nuclear and galaxy-wide 2-10 keV luminosities pro- 
vided by Young et al. (2001). Sources without listed values of 

^2-8 keV^2-8 keV na d both galactic and nuclear components 
that were undetected in the 2-8 keV band. For NGC 23, the 
nucleus was detected in the 2-8 keV band, but due to the in- 
creased background over the larger galaxy-wide area, the total 
galaxy emission was formally undetected; for this source, we 

listed the 3a lower limit on $™I keV^f-8 keV We consider 
the observed 2-10 keV emission of the LIRGs to be "AGN 
dominant" if there is both evidence for AGN activity from 

9 Feigelson et al. (2000) is available on the WWW at 
http://www.astro.psu.edu/xray/acis/memos/memoindex.html 



the optical, infrared, and/or X-ray as noted in column (10) 

in Table 3 and $™8 kev/ $ !-8 keV ^ °- 5 ( see column [9] of 
Table 3). We find that NGC 1068, NGC 1365, NGC 4194, 
and NGC 7771 satisfy these conditions. In §§5 and 6 be- 
low, we study galaxies with 2-10 keV emission dominated 
by X-ray binary populations, and therefore exclude the four 
AGN-dominant LIRGs from our further analyses. We note 
that a few of the remaining "normal" LIRGs in our sample 
harbor AGNs, which likely provide some contribution to the 
nuclear 2-10 keV emission. However, since we generally ex- 
pect that the unresolved nucleus (at scales ^ 500 pc) will also 
contain powerful star-formation regions (e.g., Scoville etal. 
2000; Alonso-Herrero et al. 2006a), which will contribute to 
the total galaxy-wide X-ray binary emission, we do not at- 
tempt to subtract any AGN component in estimating L^. 

5. RESULTS 

In this section, we explore the correlation between X-ray 
luminosity and SFR using exclusively Chandra data from the 
nearby normal galaxies from C04, LIRGs from this work, and 
LIRGs/ULIRGs from the GOALS (Iw09); hereafter, we refer 
to these three data sets collectively as the full Chandra sam- 
ple. In the analyses that follow, we have removed galaxies 
that are observed to harbor an X-ray luminous AGN, and for 
the C04 sample, we removed early-type galaxies that do not 
have accurate measures of SFR. For the C04 sample, AGNs 
played a negligible role in the total X-ray emission from their 
galaxies and therefore all 24 late-type galaxies (i.e., spirals 
and irregulars) are included here. For the LIRG sample, the 
AGNs were identified following the techniques in §4.2; the 
remaining 13 "normal" LIRGs are members of the full Chan- 
dra sample. For the Iw09 sample of 44 LIRGs/ULIRGs, we 
excluded the 15 obvious AGNs noted by Iw09. We note, how- 
ever, that Iw09 utilized X-ray hardness and/or the presence 
of a strong Fe K line at «6.4 keV to identify AGN activity 
in their LIRGs/ULIRGs, which is somewhat different from 
the methods used for our LIRG sample (see §4.2). For ex- 
ample, Arp 299, the only galaxy present in both the Iw09 
LIRGs/ULIRGs and our LIRG sample, was classified as an 
AGN by Iw09; hereafter, we adopt our "normal galaxy" clas- 
sification and physical properties for this source. For the re- 
maining 29 "normal" LIRGs/ULIRGs from the Iw09 sample, 
23 had estimates of M* and SFR from H10, and for the 6 
galaxies without H10 counterparts, we utilized the methods 
described in §2.2 to calculate M t and SFR. In total, our full 
Chandra sample consists of 66 normal galaxies (i.e., 24 from 
C04, 13 LIRGs from this study, and 29 LIRGs/ULIRGs from 
Iw09) and we use these galaxies in the statistical analyses that 
follow. 

5.1. A Chandra Perspective on the X-ray/SFR Correlation 
In Figure 4a, we show the galaxy-wide 2-10 keV luminos- 
ity ZJqx versus SFR for the full Chandra sample. Beginning 

with a simple constant L^/SFR model, logLj^/SFR = A, 
where A is a constant, we utilized the Kaplan-Meier estimator 
(e.g., Feigelson & Nelson 1986) available through the Astron- 
omy SURVival Analysis software package (ASURV Rev. 1.2; 
Isobe & Feigelson 1990; LaValley etal. 1992) to calculate the 
mean logarithmic ratio A. The Kaplan-Meier estimator ap- 
propriately handles the censored data (i.e., Lf^ upper lim- 
its) in our sample. We find A = 39.24 ± 0.06 (in logarithmic 
units of ergs s -1 [M yr -1 ] -1 ; la error) with a lcr scatter of 
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FIG. 4. — (a) Galaxy-wide 2-10 keV luminosity Lj^ versus SFR for 
the full Chandra sample: nearby normal galaxies from C04 (gray trian- 
gles), LIRGs from this study (black circles), and LiRGs/ULIRGs from the 
Iw09 GOALS study (gray squares) are plotted. We have highlighted "nor- 
mal" LIRGs in our sample that show some evidence for AGN activity (see 
column [10] of Table 3) with open circles. For reference, we have plot- 
ted the location of the MW as a cross. We have plotted the mean value of 
log Lj^ x /SFR as a dashed line and our best-fit relation assuming a power-law 

with logZ|j X = B + ClogSFR as a dotted line. The expected SFR-variable re- 
lation is plotted as a solid curve; this relation was derived using our best-fit 
relation for L^ a x = aM„ + /3SFR and the relationship between SFR and M* 
plotted in Fig. 4b. We note that this curve predicts well the general trend for 
the X-ray emission per SFR and tends to decrease with increasing SFR be- 
yond SFR es 1-10 Mq yr" 1 . (b) M k versus SFR for the full Chandra sample; 
symbols have the same meaning as they did in Fig. 4a. The dot-dashed and 
long-dashed curves highlight the best-fit relationship between and SFR 
for galaxies with SFR ^ 5 Mq yr~' and SFR ^, 5 Mq yr -1 , respectively. 
We observe that for nearby galaxies with SFR 5 Mq yr~' , M* and SFR 
are roughly linearly proportional; however, at higher SFR values, M*/SFR 
becomes smaller with increasing SFR. 

0.48 dex; this relation is plotted in Fig 4a with a dashed line 
and is tabulated in Table 4. Assuming a power-law scenario 
where logLg^ = B + Clog SFR, we utilized the expectation- 
maximization (EM) algorithm available through ASURV to 
find a best-fit relation of B = 39.46 ± 0.06 and C = 0.76 ± 0.04 
(see Fig. 4a, dotted curve and Table 4). We note that the 
use of a power-law scenario reduces the scatter by «0. 1 dex 
(«0.39 dex scatter); however, the assumption of such a highly 
non-linear relation over the entire SFR range is difficult to ex- 
plain on physical grounds. 

We note that at SFR £ 1-10 M yr" 1 , the slope of the 
X-ray/SFR relation appears to become flatter. Grimm etal. 
(2003) and Gilfanov et al. (2004a,b) argue that such a change 
in slope is expected in a scenario where the X-ray emission 



is dominated by HMXBs and a universal HMXB luminosity 
function with a normalization that scales linearly with SFR is 
applied. In this scenario, at SFR ;$ 1-10 M Q yr" 1 , the rela- 
tion is expected to be steeper than linear (Z^x x SFR 1 7 ), and 
at SFR £ 1-10 M yr" 1 , the relation is expected to become 
linear (Z^x x SFR). Dividing our data into low and high SFR 
regimes around SFR = 5 M Q yr" 1 and performing power-law 
fits, we find 



, g ai _ I 10( 39 - 57±011 )SFR (0 - 94±015) SFR < 0.4 M yr" 1 

^HX - i 10 (39.49±0.21) SFR (0.74±0.12) spR > g 4 M& yJ .-l ^ 



where the cut-off SFR « 0.4 M yr 1 was determined by set- 
ting the upper and lower relations for equal to each other 
and solving for SFR. This two-regime relation can be used to 
obtain a reasonable (with a lcr scatter of «0.4 dex; see Ta- 
ble 4) estimate of L^ x based on SFR alone; however, as we 
will discuss in § 6.2 below, some LIRGs and ULIRGs may 
fall «0.5-l dex below this relation due to significant intrinsic 
attenuation of the 2-10 keV emission. We note that the mea- 
sured slopes of these relations are inconsistent with those pre- 
dicted by Grimm etal. (2003) and Gilfanov etal. (2004a,b). 
We speculate that this is due to the presence of LMXBs that 
provide a non-negligible contribution to in the low-SFR 
regime. In the next section, we explore the use of a physically- 
motivated scaling of L™ with M t and SFR that incorporates 
emission from both HMXB and LMXB populations. 

5.2. Correlating Lf^ with M t and SFR 

From Figure 4a, we notice that the mean L^/SFR ratio 
(dashed line) for the total Chandra sample generally lies be- 
low the data for galaxies with SFR ^ 5 M Q yr" 1 (i.e., from 

C04), primarily due to the lower values of L^ d x /SFR for the 
many LIRGs/ULIRGs in the sample. This trend has been 
noted by PR07, who suggest that typical nearby galaxies 
(SFR ;S 5 Mq yr" 1 ) will have fractionally important contri- 
butions from both HMXBs and LMXBs, making them more 
X-ray luminous per unit SFR than "star-formation active" 
galaxies like LIRGs/ULIRGs, which are likely to have neg- 
ligible contributions from LMXBs and will be dominated by 
HMXBs. As recognized by C04, the total 2-10 keV luminos- 
ity for nearby galaxies can reasonably be quantified as the 
sum of contributions from LMXBs and HMXBs, which to 
first order are expected to scale linearly with M t and SFR, 
respectivelyCB The C04 relation can be expressed using the 
following equation: 



r gal 



Lg^(LMXB) + Lg"x( HMXBs ) = aM * + /3SFR, (3) 

where a and /3 are scaling constants (see below). From equa- 
tion 3, we expect that galaxies with larger SFR/M* will have 
more significant contributions from HMXBs, whereas galax- 
ies with smaller SFR/M* will have larger contributions from 
LMXBs. We would therefore expect that should scale 
more closely with SFR/M+ than SFR alone. 

In Figure 4b, we show versus SFR for the full Chan- 
dra sample. It is apparent that there is a strong corre- 
lation between M t and SFR (Spearman's p probability of 
>99.99%). Using simple polynomial fitting techniques, 

10 However, see Grimm etal. (2003), Gilfanov etal. (2004a), and Gil- 
fanov (2004) for different scalings depending on M„ and SFR regimes. 
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FIG. 5. — (a) Logarithm of the galaxy-wide 2-10 keV luminosity per SFR 
(Zjj^/SFR) versus SFR/M* for the full Chandra sample; symbols are the 
same as they were in Fig. 4a. The PR07 relations for all galaxies and HMXBs 
are shown as dashed and dotted lines, respectively. The C04 relation and our 
best-fit relation are shown as dashed-dot and solid curves, respectively. The 
clear dependence of Lg^/SFR on SFR/M* convincingly illustrates that SFR 
plus provides a superior estimate of than a simple L^/SFR correla- 
tion, (b) Residuals to the best-fit relation = aM 4 +/3SFR. Horizontal dot- 
ted lines show the \a scatter for this relation. We note that for the most star- 
formation active galaxies, there are several sources that are X-ray underlumi- 
nous; we have annotated the locations of NGC 4418 and IRAS F10173+0828, 
the most extreme outliers of the LIRGs and ULIRGs, respectively. For such 
sources, it is plausible that significant extinction (Ay ^ 50 mag), which has 
been observed in powerful ULIRGs (see § 6.2), will play a role in obscur- 
ing the 2-10 keV emission. For reference, we have plotted the expected 
attenuation of the 2-10 keV emission for optical extinction in the range of 
A v = 0-500 mag (corresponding to N K sa [0-1.2] X 10 24 cm" 2 ). 



we find for galaxies with SFR <, 5 M Q yr" 1 , M t and SFR 
are nearly linearly correlated (M* oc SFR 1 ' 10 ); however, 
for galaxies with SFR ^ 5M Q yr _1 , the relationship be- 
comes shallower (M* oc SFR ' 7 ). From this analysis, we 
find that the quantity SFR/M* remains roughly constant for 
SFR w 0.01-1 Mq yr -1 , and to first order, this plausibly sug- 
gests that the HMXB and LMXB fractional contributions 
would remain roughly constant as well. Indeed, as noted in 
§ 5.1 above, we found that Z|£ oc SFR (0 ' 94±015) over this 
range of SFR. At SFR > 5 M & yr -1 , the average SFR/M* 
increases from a constant value, which would qualitatively al- 
low for HMXBs to make an increasingly large fractional con- 



tribution to the total X-ray luminosity while the Lj[ X /SFR ra- 
tio decreases due to the reduced fractional contribution from 
LMXBs; in § 5.1 above, we found that Lg£ oc SFR (0 ' 74±012) 



in this regime. 

The best estimates to date of the quantities a and j3 come 
from the C04 analysis, which does not include Chandra con- 
straints on galaxies with SFR £ 5 M Q yr" 1 (i.e., LIRGs and 
ULIRGs), a regime where the SFR/M* is expected to be large 
and HMXBs are likely to dominate. Here, we utilize the full 
Chandra sample to improve the calibration of a and (3 for data 
covering the broad range of SFR « 0.01^-00 M Q yr" 1 . 

Considering the above discussion, it is useful to examine 
the relationship of L^/SFR as a function of SFR/M*. Fol- 
lowing equation 3, the quantity L^ X /SFR can be written as 
follows: 

Lgx/SFR = (aM* + /3SFR) /SFR = aCSFR/M*)" 1 + /3. (4) 

From equation 4, we see that Lg^/SFR oc (SFR/M*)" 1 , and 

in Figure 5a, we plot Lj^ x /SFR versus SFR/M* for the full 
Chandra sample. We see that the IMF and X-ray bandpass 
adjusted C04 relation (dot-dashed curve in Fig 5a), which 
was determined solely from the C04 data (triangles), pre- 
dicts well the basic observed trend; however, it mildly un- 
derpredicts Lfx/SFR for LIRGs and ULIRGs. Using the full 
Chandra sample, we fit the data using the EM algorithm and 
the model provided in equation 4. We find best-fit values of 
a = (9.05 ± 0.37) x 10 28 ergs s" 1 Mq 1 and /? = (1.62 ± 0.22) x 
10 39 ergs s" 1 (M Q yr" 1 )" 1 with a scatter of 0.34 dex, which is 
smaller than that produced by fitting for SFR alone (see §5.1; 
Table 4). 

Our best-fit relation is plotted as a solid curve in Figure 5a. 
In Figure 4a, we combine this relation with the M*-SFR re- 
lations provided in § 5.1 and Figure 4b to show the expected 
Lj^-SFR relation (see solid curve in Fig. 4a); in general, the 
combined relations appear to fit the data better than a sim- 
ple mean logLj^/SFR ratio and a less physically meaningful 

power-law relation (i.e., logL^ x = B + CSFR). 

After accounting for differences in IMF and X-ray band- 
pass, we find that a/a(C04) = 0.76 and /3//3(C04) = 2.18, 
suggesting that for a given SFR/M+, HMXBs provide «2.9 
times larger L x ' contributions over LMXBs than previously 
reported. We further note that our best-fit relation has reduced 
the error bars for a and j3 by factors of «3.8 and «2.1, re- 
spectively. We find that for SFR/M* > 5.6 x 10" 11 yr -1 (cor- 
responding roughly to SFR > 2 M yr" 1 or L IR > 2 x 1O 1O L 
assuming the M*-SFR correlation in Fig. 4b and 7 = 1 in equa- 
tion 2), HMXBs will dominate the galaxy-wide 2-10 keV. 
This transition SFR/M* is a factor of «2.9 times lower than 
previous estimates. 

We note that the MW (plotted in Fig. 5 as a cross) ap- 
pears to be located below our best-fit correlation line by a 
factor of «0.37 dex, which is just outside the ler scatter 
in the relation. Our best-fit relation predicts that a galaxy 
with M± and SFR equal to those of the MW, will have 
Lg x (LMXB)/Lg x (HMXB) w 1 .4, which is much less than the 

Lg x (LMXB)/Lg x (HMXB) - 10 estimated by Grimm etal. 
(2002) using the measured luminosity functions of LMXBs 
and HMXBs. We speculate that these differences are due to 
either (1) the large uncertainties in measuring M + and SFR 
for the MW, as a result of our location within it, and/or (2) 
the MW being a true outlier to the general correlation. It has 
been noted by Hammer et al. (2007), that the MW does not 
appear to follow other correlations along physical parameter 
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TABLE 4 
Best-Fit X-ray/SFR Relations 



Model Description 

logLg 1 ' =A+logSFR 

logZ|| = .B+ Clog SFR 

log£g|=B+ClogSFR(SFR <O.4M yr') 
logLfif =B + ClogSFR (SFR > 0.4 M yr 1 ) 
L^ = aM± + /3SFR 



Best-Fit Parameters 



A = 39.24 ±0.06 
B = 39.46 ± 0.06 and C = 0.76 ± 0.04 
B = 39.57±0.11 and C = 0.94 ±0.15 
B = 39.49 ±0.21 and C = 0.74 ±0.12 
a = (9.05 ± 0.37) X 10 28 ergs s" 1 M Q ' 
g = (1.62 ± 0.22) X 10 39 ergs s~' (M yr" 1 ) 



lo- Scatter (dex) 



0.48 
0.39 
0.41 
0.36 

0.34 



NOTE. — Summary of best-fit relations between SFR, and for the models described in § 5. Fits have been performed following the methods described 
in §§5.1 and 5.2; all quoted errors are la. We note that the inclusion of both SFR and M* in the final model provides the least statistical scatter of all the models. 



planes (e.g., planes drawn by disk rotation velocity, /^-band 
luminosity, disk angular momentum, and mean [Fe/H] abun- 
dances of stars in the outskirts of the disk) drawn using a sam- 
ple of galaxies with similar M*. Hammer etal. conclude that 
the MW is likely to have a different star-formation history 
and disk construction than other similar galaxies; such condi- 
tions could plausibly affect the scaling of LMXB and HMXB 
power output with and SFR, respectively. 

6. DISCUSSION 

In Figure 5b, we show the residuals to our best-fit relation 
(i.e., L x /[aM+ + (3SFR]). As noted in § 5.2, the scatter over 
the entire range of SFR/M* is «0.34 dex. There are several 
factors that are likely to contribute to the scatter in this rela- 
tion including: (1) variations in the LMXB populations due 
to a diversity of stellar ages and star-formation histories (e.g., 
Fragos et al. 2008); (2) unrelated low-level AGN activity that 
affect estimates of the X-ray and/or infrared luminosities, two 
quantities that are important to this analysis (e.g., Imanishi 
et al. 2007, 2010; Iw09); (3) strong statistical fluctuations due 
to small numbers of luminous X-ray binaries that dominate 
the total X-ray luminosity, particularly for galaxies with small 
SFR and/or (e.g., Grimm et al. 2003; Gilfanov et al. 2004a; 
see § 6.1 below); and (4) significant extinction that can affect 
estimates of stellar mass, and in the most extreme cases, may 
obscure even 2-10 keV emission (see § 6.2 below). In §§ 6.1 
and 6.2 below, we explore in more detail points (3) and (4), 
respectively. 

6.1. Statistical Considerations at Low SFR and M t 

As discussed in § 5.1, the qualitative change in slope of the 
X-ray/SFR correlation has previously been noted by Grimm 
et al. (2003) and Gilfanov et al. (2004a,b) who predict such 

behavior for an HMXB -dominant Lf^. These authors suggest 
that since the HMXB luminosity function is a relatively shal- 
low power-law with dN/dLxp oc SFR L^p (with 7 « 1.6) and 
a maximum cut-off luminosity of L™p~° ff ~ 2 x 10 4() ergs s" 1 
(see, e.g., Grimm etal. 2003; Swartz etal. 2004), then in the 
low-SFR regime, the most luminous HMXB will have a most 
probable luminosity L™p x oc 

SFR i/(7 ; which is less than 

L™p" off . In this regime, the integrated X-ray luminosity Lf^ 

is therefore expected to scale with SFR I//(7_I) (or SFR 1 7 in 
the case of 7 = 1.6). Such an effect is predicted to produce 
a non-linear SFR dependence for the HMXB contribution to 
Z|jx for SFR 1-10 M Q yr" 1 . Similar logic can be applied 

to the LMXB contribution to using the more complex 
luminosity function for LMXBs, which is expected to have 
a linear normalization that scales with M*. However, the 
cut-off luminosity for the LMXB population is likely to be 



much lower (^10 39 ergs s ') than that of HMXBs and there- 
fore Lf[x(LMXB) is expected to be linear for galaxies with 
£ 5-20 x 10 9 M© and very close to linear for smaller val- 
ues of (e.g., Gilfanov etal. 2004a; Gilfanov 2004). 
In § 5.1, we noted that in the low-SFR regime (SFR w 

0.01-1 M© yr" 1 ), Lf^/SFR and SFR/M* were roughly con- 
stant; over this range these quantities have median values 
of «5 x 10 39 ergs s _1 (M Q yr -1 ) -1 and 8 x 10" 11 yr -1 , re- 
spectively. From our best-fit relation for equation 3, we es- 
timate L^LMXByLg^HMXB) w 0.74 for these galaxies, 
suggesting that both LMXBs and HMXBs are estimated to 
provide significant contributions to in this regime. How- 
ever, if indeed Lg^HMXB) oc SFR 1 7 and Lg£(LMXB) <* 
M* oc SFR are accurate descriptions of the respective HMXB 
and LMXB emission in this regime, then we would have ex- 
pected that Lg£ oc SFR 1 0-1 7 



In § 5.1, we found that oc 



SFR (094±0 15) . This suggests that either (1) Lf^(HMXB) 

scales with SFR 1 7 and provides a small contribution to Lf^ or 

(2) Lg^HMXB) <* SFR - Future X-ray and multiwavelength 
studies of large numbers of normal galaxies in this regime 
may provide the necessary statistical basis to allow discrim- 
ination between these two scenarios; future X-ray missions 
such as IXO and WF XT will enable efficient detection of many 
galaxies in the low-SFR regimeO Otherwise, a combina- 
tion of X-ray and optical/near-IR observatories with highly- 
precise astrometry would be needed to allow for the detection 
of X-ray binaries and the reliable identification and classifica- 
tion of their stellar counterparts in external galaxies. This is 
currently broadly possible with the combination of Chandra 
and HST data for a small number of nearby galaxies. Future 
missions such as Generation-X and JWST will undoubtedly 
expand this science to more representative galaxies that cover 
the entire low-SFR regimePl 

6.2. A Deficit of X-ray Emission in the Most Actively 
Star-Forming Galaxies 

From Figure 5b, it is evident that for several of 
the most star-formation active LIRGs/ULIRGs (i.e., 
SFR/M* > 10" 9 yr" 1 ) in the total Chandra sam- 
ple there is a departure from the best-fit relation. 
For example, IRAS 07251-0248, IRAS F15250+3608, 
IRAS F10173+0828, IRAS 21101+5810, and NGC 4418 

11 More information regarding IX O and WFXT can 

be found at |http : //ixo . gsf c . nasa . gov"7] and 

|http : //wf xt . pha . jhu . edu/| respectively. 

12 More information regarding Generation-X and JWST can be 
found at http://www.cfa.harvard.edu/hea/genx/ and 
http : / / www . jwst . nasa . gov/, respectively. 
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Fig. 6. — Residuals to the best-fit relation = qM* +/3SFR versus the IRAS 60-to-100/im flux density ratio ^60/^1 00 («) and the logarithm of the infrared- 
to-UV luminosity ratio logLiR/Luv (b). Symbols and lines have the same meaning as they did in Figure 5b. A Spearman's p test reveals that the residuals are 
negatively correlated with F6o/^loo and logLiR/Luv at the Ri99.8% and 99.9% confidence levels, respectively. The quantity /^o/fioo is a proxy for temperature 
and infrared luminosity density (star-formation density) and logLm/Luv is a direct measure of UV obscuration. These relations suggest that the most X-ray 
deficient LIRGs and ULIRGs, which have relatively warm dust temperatures and large UV obscuration due to dust, contain dense buried star-forming regions 
with attenuated 2-10 keV emission. 



all have residuals to our best fit that fall below the 2er 
scatter, when only ssl-2 galaxies are expected statistically 
over the entire SFR/M* range. If we exclude galaxies with 
SFR/M* £ 10~ 9 yr" 1 , this significantly reduces the scatter 
to «0.26 dex (i.e., by «20%). Studies of LIRGs/ULIRGs 
have found that in these objects the powerful star-formation 
activity occurs in compact regions that are surrounded by 
large columns of gas and dust (e.g., Soifer etal. 2000; 
Alonso-Herrero etal. 2006b; Siebenmorgen etal. 2008). The 
obscuration to these star-forming regions can be very signifi- 
cant (A v >, 50 mag; e.g., Genzel etal. 1998; Siebenmorgen & 
Kriigel 2007) and can attenuate emission over the 2-10 keV 
bandpass. 

To estimate the effect that obscuration can have on 
the observed 2-10 keV power for ULIRG levels of ex- 
tinction, we first converted Ay to A^h following the rela- 
tion N H « 2.3 x 10 21 A V cm" 2 from Giiver & Ozel (2009), 
which was calibrated using measurements of various columns 
through the Milky Way. We note that this calibration has been 
calibrated only over the range of A v « 0.2-30 mag and is 
subject to large uncertainties when applied to other galaxies, 
which may have different chemical compositions in their in- 
terstellar mediums. A range of A v = 0-500 therefore implies 
column densities in the range A^h ~ 0-1 .2 x 10 24 cm" 2 . Using 
PIMMS and an assumed power-law X-ray SED with V = 1 .8, 
we find that the emergent 2-10 keV luminosity will be atten- 
uated by factors of «2, 3, and 15 for A v w 50, 100, and 500, 
respectively (see vertical bar in Fig. 5b); such visual extinc- 
tions are likely to be present in at least some LIRGs/ULIRGs 
(e.g., Genzel etal. 1998 estimate that some ULIRGs may 
even reach Ay w 1000). We note that there are only a few 
LIRGs/ULIRGs in our sample with the high-quality optical 
to far-IR data required to estimate A v accurately. For exam- 
ple, the well-studied ULIRG Arp 220 has a 2-10 keV lumi- 
nosity that is a factor of ?»3 times lower than that expected 
from our best-fit relation for equation 4. Detailed modelling 
of the infrared-to-optical spectrum suggests a global value of 
A v w 50-150 mag (e.g., Siebenmorgen & Kriigel 2007), im- 



plying Nh w 1.2-3.5 xlO 23 cm and therefore an expected 
factor of K.7.-A attenuation of the intrinsic 2-10 keV luminos- 
ity, in good agreement with the deficit that we observe. 

To assess whether there is broader evidence that the X-ray 
underluminous LIRGs/ULIRGs in the full Chandra sample 
may be subject to significant 2-10 keV attenuation, we uti- 
lized two informative quantities: (1) the IRAS 60-to-lOO/im 
flux ratio (Feo/^iooX which has been shown to correlate well 
with mean dust temperature and infrared luminosity surface 
density (e.g., Chanial et al. 2007); and (2) the IR-to-UV lu- 
minosity ratio (Lir/Ltjv), which provides a direct measure 
of UV obscuration from star-forming regions. If the X-ray 
emission from star-forming regions is buried under significant 
absorbing columns, we would expect the star-formation to 
be highly concentrated and not dispersed through the galaxy 
(i.e., large Feo/^ioo) and heavily obscured in the UV (i.e., 
large Lir/Luv)- In Figure 6a and 6b, we plot the residuals to 
our best-fit relation (Lf^/faM* + /3SFR]) versus /^o/^ioo and 
Lir/L\jv, respectively. It is evident that the LIRGs/ULIRGs 
with the largest deficit of X-ray emission with respect to our 
best-fit relation are those with large values of F$o /F100 and 
Lir/Luv', a Spearman's p test reveals that the residuals are 
correlated with each quantity at the 99.8% and 99.9% signif- 
icance levels, respectively. We therefore conclude that ob- 
scuration is at least partially responsible for the deficit of ob- 
served 2-10 keV emission at the largest SFR/M*. 

7. CONCLUSIONS 

In this paper, we present new results from Chandra obser- 
vations of a complete sample of 17 LIRGs with D < 60 Mpc 
and Galactic column densities A^h ~ 5 x 10 20 cm" 2 . Our 
Chandra observations reveal a variety of X-ray morphologies 
for the LIRGs, with the majority of their corresponding X-ray 
spectra being well-described by a hot-gas component (with 
kT 0.8 keV) plus a power-law (r w 1-3), which is expected 
to originate from X-ray binaries, and when present, an AGN. 
We find that four out of the 17 LIRGs both show evidence for 
AGN activity in the optical, infrared, and/or X-ray band and 
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contain a nuclear component that dominates the integrated 
2-10 keV luminosity. We study the galaxy-wide 2-10 keV 
luminosity Lf^ and its correlation with SFR and M t using 
the 13 non-AGN-dominant LIRGs combined with Chandra 
studies of less star-formation active normal galaxies (from 
C04) and more powerful LIRGs and ULIRGs (from Iw09), 
collectively referred to as the total Chandra sample. Our key 
results can be summarized as follows: 

(1) — For the galaxies that make up our total Chandra sample, 
we find that SFR and M* are strongly correlated, such that 
M* oc SFR 1 10 for SFR w 0.01-5 M Q yr" 1 and M* oc SFR 27 
for SFR > 5 M Q yr" 1 (see Fig. 4b). 

(2) — We find that is strongly correlated with SFR (the 
X-ray/SFR correlation), as has been found in previous stud- 
ies (see Fig. 4b). Similar to point (1), the X-ray/SFR cor- 
relation slope becomes shallower for SFR £ 1-10 M Q 

yr" 1 . We find that Lg£ = (39.57 ± 0.1 l)SFR (a94±(u5) for 
SFR < 0.4M o yr-' and Lg£ = (39.49 ± 0.21)SFR (0 ' 74±012) 
for SFR £ 0.4 M Q yr" 1 provides a reasonable (with w0.4 dex 
scatter; see Table 4) estimate of the L^y. The change in slope 
is likely to be related to the changing slope of the correla- 
tion between SFR and M*, since we expect the total galaxy- 
wide 2-10 keV luminosity to be the sum of contributions from 
HMXBs and LMXBs, which to first order are expected to 
scale with SFR and M*, respectively. 



(3) — Under the assumption that L^ x can be written as L^ x = 
Lg'x(LMXB) + Lf^(HMXB) = aM, + /3SFR, we find best- 
fitting parameters of a = (9.05 ±0.37) x 10 28 ergs s" 1 Mq 
and p = (1.62 ±0.22) x 10 39 ergs s" 1 (M Q yr" 1 )" 1 . These 
values suggest that HMXBs dominate the X-ray emission for 
SFR/M* > 5.9 x 10" 11 yr" 1 , a factor of «2. 9 times lower than 
previous estimates. 

(4) — For the most actively star-forming LIRGs and ULIRGs 
(i.e., those with SFR/M* £ 10~ 9 yr" 1 ), we find an excess of 

objects with lower than expected. These sources have 
infrared colors and UV to total-IR luminosities indicative of 
compact star-forming regions with extreme extinction. We 
argue that such star-forming regions are likely to have attenu- 
ated 2-10 keV emission. 

We thank the referee for their helpful comments that have 
improved the manuscript. We thank Kazushi Iwasawa for 
generously sharing data and Lee Armus and Joseph Maz- 
zarella for helpful discussions. We gratefully acknowl- 
edge financial support from the Einstein Fellowship Pro- 
gram (B.D.L.), the Royal Society (D.M.A.), the Leverhulme 
Trust (D.M.A.), Chandra X-ray Center grant G09-0134A/B 
(W.N.B.,F.E.B.), and NASA ADP grant NNX10AC99G 
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